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Dynamical few-body encounters in the dense cores of young massive star clusters are responsible for the loss of a significant fraction 
of their massive stellar content. Some of the escaping (runaway) stars move through the ambient medium supersonically and can be 
revealed via detection of their bow shocks (visible in the infrared, optical or radio). In this paper, which is the second of a series 
of papers devoted to the search for OB stars running away from young (^ several Myr) Galactic clusters and OB associations, we 
, present the results of the search for bow shocks around the star-forming region NGC 6357. Using the archival data of the Midcourse 

. Space Experiment (MSX) satellite and the Spitzer Space Telescope, and the preliminary data release of the Wide-Field Infrared Survey 

Explorer (WISE), we discovered seven bow shocks, whose geometry is consistent with the possibility that they are generated by stars 
r"| \ expelled from the young (~ 1-2 Myr) star clusters, Pismis24 and AH03 J1725-34.4, associated with NGC 6357. Two of the seven 

O ^. bow shocks are driven by the already known OB stars, HD 319881 and [N78] 34. Follow-up spectroscopy of three other bow-shock- 

I . producing stars showed that they are massive (0-type) stars as well, while the 2MASS photometry of the remaining two stars suggests 

O ' that they could be BO V stars, provided that both are located at the same distance as NGC 6357. Detection of numerous massive stars 

^ I ejected from the very young clusters is consistent with the theoretical expectation that star clusters can eff'ectively lose massive stars 

C/) , at the very beginning of their dynamical evolution (long before the second mechanism for production of runaway stars, based on a 

. supernova explosion in a massive tight binary system, begins to operate) and lends strong support to the idea that probably all field 

OB stars have been dynamically ejected from their birth clusters. A by-product of our search for bow shocks around NGC 6357 is the 
_ detection of three circular shells typical of luminous blue variable and late WN-type Wolf-Rayet stars. 

^ ' Key words, stars: kinematics and dynamics - stars: massive - open clusters and associations: general - open clusters and associations: 

, individual: Pismis 24 - open clusters and associations: individual: AH03 J1725-34.4 - ISM: individual objects: NGC 6357 

■ 1 ■ Introduction Hohle [2QTT1) , or indirectly, through the detection of bow shocks 

^ , , , . , . , , r around them (Gvaramadze & Bomans I2008bt Gvaramadze et 

O . Close few-body dynamical encounters in the dense cores of [2oToa Gvaramadze, Kroupa & Pflamm-Altenburg l20T0dl 

— young massive star clusters are responsible for the loss of a sig- Gvaramadze, Pflamm-Altenburg & KroupaEOlIi Gvaramadze 

nificant fraction of OBstars at the early stagesof clusterevo- al. l20TTB . The latter possibility is especially helpful for those 

ution (Poveda et al^^l Aarseth&Hills [TMl Gies [1983 runaways whose proper motions are still not available or are 

Leonard & Duncan 1990; KroupaOT Pflamm-Altenburg & j^easured with a low significance. The geometry of detected 

Kroupa 2QQ6; Moeckel & Bate 2010). The high central densi- ^^^^ -^^^^ ^j^^^^j^^ ^j^jl^^ 

H ties in young clusters (the necessary condition for the produc- ^^^-^^ (^^^^1, determine possible parent clusters for 

■ - tion of sta'-s) be either primordial (e.g O^e & ^j^^ bow-shock-producing field stars (Gvaramadze & Bomans 

PringleHMl Murray & Lin J996; Clarke & BonnellEOOi or Gvaramadze et al l20T051d.l20m. 
caused by dynamical mass segregation (e.g. Portegies Zwart et 

al. 119991 Giirkan et al. 120041 Allison et al. 12009b . In both cases, The present paper is the second of a series of papers devoted 

the clusters start to eject stars long before the most massive clus- to the search for OB stars running away from young (^ sev- 

ter members explode in supernovae, i.e. long before the binary- eral Myr) Galactic clusters and OB associations. In Gvaramadze 

supernova ejection mechanism (Blaauw 1961) begins to operate. & Bomans ( 2008b; hereafter Paper I), we reported the detec- 

Moreover, runaway stars could leave their parent clusters already tion of three bow shocks produced by OB stars running away 

during the cluster formation process if a core of massive proto- from the star cluster NGC 6611. Now we report the discovery 

stars forms before bulk gas expulsion from the embedded cluster, of seven bow shocks within ~ 1?5 from the star-forming region 

The runaway OB stars can be revealed either directly, via NGC 6357 and its associated (massive) star clusters Pismis 24 

measurement of their proper motions and/or radial velocities and AH03 J1725-34.4. The geometry of all seven bow shocks 

(e.g. Moffat et al. |1998t Mdzinarishvili & Chargeishvili 120051 suggests that they are driven by stars expelled from NGC 6357. 

Massey et al. 120051 Evans et al. 120101 TetzlafF, Neuhauser & In Sect.[2lwe review the relevant data for the star clusters associ- 
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ated with NGC 6357. Sect.[3]presents the results of the search for 
bow shocks around NGC 6357. The bow- shock-producing stars 
are identified and discussed in Sect.lH In Sect.[5]we present the 
discovery of three circular shells around NGC 6357. Sect.[6]deals 
with questions related to the content of the paper. We summarize 
in Sect.H 



2. Star-forming region NGC 6357 and its associated 
star clusters 

NGC 6357 is an extended (~ 40' x 60' or ~ 20 pc x 30 pc at 
a distance of 1.7 kpc; see below) Hii region of ongoing mas- 
sive star-formation in the Sagittarius arm. At radio wavelengths 
NGC 6357 is dominated by the two components G 353.2-^0.9 
and G353.1-h0.6 (Schraml & Mezger[M2l FelH et al. [T990b . 
Two open star clusters, Pismis 24 and AH03 J1725-34.4, are as- 
sociated with these components, respectively. The first cluster 
was extensively studied during the last years, while only sparse 
information exists on the second one. Below we review the ob- 
servational data on both clusters in turn. 



2.1. Pismis 24 

Pismis 24 was recognized as an open star cluster by Pismis 
(119591 ). It is believed that the massive stars of this cluster are the 
main ionizing source responsible for the origin of the H ii region 
NGC 6357 (Lortet, Testor & Niemela[T98l Bohigas et al. 120041 
Cappa et al. 1201 11) . Line observations of a molecular cloud asso- 
ciated with NGC 6357 showed that most of the molecular emis- 
sion arises from the regions behind or to the north of Pismis 24, 
which indicates that the cluster is immersed in a blister-like H ii 
region viewed face-on (Massi, Brand & Felli |1997b . 

The cluster contains two very massive 03.5-type stars 
(Massey, DeGioia-Eastwood & Waterhouse 120011 Walborn et 
al. 120021) , one of which, Pismis 24-1, was believed to be one of 
the most massive known stars in the Galaxy, with an inferred 
mass larger than 200 M© (Walborn et al. 120021) . Subsequent 
Hubble Space Telescope observations resolved Pismis 24-1 into 
two visual components, while follow-up spectroscopy showed 
that one of the components (Pismis 24-lSW) is an 04 III(f+) star 
(of mass of ~ 100 M©) and the second one (Pismis 24-1 NE) is 
a very massive (~ 100 M©) short-period spectroscopic binary 
(Maiz Apellaniz et al. .2007 ). 

Because of the high visual extinction towards Pismis 24 (^ 
5-15 mag; Wang et al. 2007), the stellar content of the cluster 
was poorly studied until recently, so that only the most lumi- 
nous and less reddened cluster members (about two dozen alto- 
gether) were identified, either by means of photometry or spec- 
troscopy (Moff'at & Vogt[T973t Neckel [T978l[T984l Lortet et al. 
119841 Massey et al. 120011) . One of these stars is a binary system 
(HD 157504) composed of a WC7 Wolf-Rayet star (WR 93) and 
an 07-9 star (van der Hucht [200T]) . 

Thanks to the ability of X-ray observations to pen- 
etrate heavy absorption, the situation was improved with 
Chandra/ ACIS observations of NGC 6357, which allowed the 
detection of even the low-mass pre-main- sequence population 
of Pismis 24 (with masses extending down to ~ 0.3 M©), thereby 
increasing the number of the cluster members by a factor of ~ 50 
and bringing it to ^ 10000 (Wang et al. [20071) . Most of these 
stars are concentrated in a compact core of angular radius of ^ 2 
arcmin (or ^1.0 pc), which is centred on the two most massive 



stars in the cluster, Pismis 24-1 and Pismis 24-1^1, while the re- 
maining stars are spread with an approximately exponential dis- 
tribution over a halo of radius of ^ 11 arcmin (^ 5.4 pc). The 
Chandra observation also detected two dozens of X-ray sources, 
whose luminosities and colours suggest that they could be OB 
stars. If spectroscopic follow-ups of these stars will prove that 
they are indeed massive, then the OB star content of Pismis 24 
will be doubled. Assuming that Chandra detected all stars with 
mass > 0.3 M© and using the Kroupa (12001b IMF (which is a 
two-part power-law IMF with a slope ai = 1 .3 for stellar masses 
between 0.08 and 0.5 M©, and a Salpeter slope a2 = 2.3 for 
more massive stars), one finds a mass of the cluster of ~ 10"^ M© 
and the expected number of OB stars in the cluster of ~ 100. 
The latter figure exceeds the number of confirmed and candi- 
date OB stars (i.e. 44 stars; Wang et al. 120071) by a factor of 2, 
from which one can conclude that either the IMF of the cluster's 
members is steeper than the Salpeter one or most OB stars were 
already ej ected f rom the cluster (cf. Paper I; Pflamm-Altenburg 
& Kroupa l20Q6l Moeckel & Batel20T0n In the latter case, one 
can expect to find numerous OB stars all around the cluster (see 
Sect.O. 

The simultaneous presence in Pismis 24 of several very 
young (~ 1 Myr) and very massive (^ 90 - 100 M©) stars (Maiz 
Apellaniz et al. l2007l) and the more evolved (^ 2 — 2.5 Myr 
old) WC7 star WR93 (Massey et al. l200T1) leads to an age dis- 
crepancy. To overcome this discrepancy, one can assume that 
the star-formation process in the cluster was not coeval or that 
WR 93 is projected by chance near the line of sight to the cluster 
(note that WR 93 is located only 4 arcmin away from the clus- 
ter centre, i.e. well within the cluster halo). In the latter case it 
is plausible that the parent cluster of WR 93 is associated with 
the same molecular cloud as Pismis 24 but is totally obscured by 
the dusty material of the cloud, while the star itself is observ- 
able only because it is a runaway, which already escaped from 
the densest part of the cloud (cf. Gvaramadze et al. 2010a'). Note 
that the coexistence of evolved and very young massive stars in 
the same star-forming region is not a rare occurrence. For ex- 
ample, two evolved (~ 4 Myr old) massive stars are projected 
against the young (^ 1 Myr) very massive cluster NGC 3603, 
one of which, the blue supergiant star Sher25, is separated from 
the cluster core by less than 20 arcsec (or 1 pc in projection; e.g. 
Melena et al. I2008[ Rochau et al. 120101) . Another possible ex- 
planation for the age discrepancy is that the most massive stars 
in the cluster are blue stragglers, formed because of merging of 
less massive stars in the course of close dynamical encounters in 
the cluster's core (cf. Gvaramadze & Bomans ?; Paper I). 

Numerous distance estimates for Pismis 24 range from ^1.7 
to 2.6 kpc (e .g. Moff'at & VogtWT^! Neckel lT978llT984l Lortet 
et al. I1984[ Massey et al. 2001j , putting the cluster in the 
Sagittarius spiral arm. To constrain the distance to Pismis 24, we 
used the optical (Massey et al. 1200 II) and near-infrared (2MASS; 
Skrutskie et al. 120061) photometry of all (six) known dwarf O 
stars in the cluster (see Table 2 in Massey et al. 120011) and the 
UBVJ HK syn thetic photometry of Galactic O stars by Martins 
& Plez (120061) . Assuming that the standard total-to- selective ab- 

^ Note that the coordinates of Pismis 24 given in the SIMBAD 
database are significantly (~ 16 arcmin) off the position of these stars 
and rather correspond to the second star cluster, AH03 J1725-34.4, as- 
sociated with NGC 6357. The latter cluster is confused with Pismis 24 
in the WEB DA database (http://www.univie.ac.at/webda/) as well. 

^ The discrepancy between the "observed" and expected number of 
OB stars in Pismis 24 would be less severe if most massive stars residing 
in the cluster are binaries with mass ratio close to unity (cf. Clarke & 
Prinde [T992l ). 
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sorption ratio Ry = Av/E(B -V) = 3. 1 is valid for Pismis 24 and 
using the optical photometry, one finds a true distance modulus 
of ^ 12.0 mag and a distance of 2.5 kpc, which agree well with 
those derived by Massey et al. (12001b on the basis of a diff'erent 
calibration of stellar parameters. On the other hand, adopting the 
extinction law from Rieke & Rebofsky (11985b , so that 

Ak,=0.66E(J-K,), (1) 

and using the 2MASS photometry, one finds a true distance mod- 
ulus of ^ 11.15 mag and a distance of 1.7 kpc. The discrep- 
ancy between the two distance estimates could be understood if 
the reddening towards Pismis 24 is anomalous, i.e. Ry > 3.1. 
One can show that Ry - 3.6 would be enough to reconcile the 
distance estimates (cf . with Neckel & Chini 119811 who found 
Ry = 3.8 for NGC 6357). The distance of 1.7 kpc could also be 
derived if one estimates the visual extinction Ay using the rela- 
tionship 

Ay = 139E(V-J), (2) 

which, according to Tapia et al. (I1988K is valid across the 
Galaxy, irrespective of high Ry values (see also Persi & Tapia 
120081) . This distance is consistent with that to the Hii region 
NGC 6334 1.7 kpc; Persi & Tapia|2008l), located on the sky 
at ~ 2° to the southwest from NGC 6357. It is beUeved that both 
these H ii regions form a single high-mass star-forming complex 
(Russeil et al. 120 101 and references therein). In the following we 
adopt a distance of 1.7 kpc for Pismis 24 (and the star-forming 
region NGC 6357 as a whole) so that 1° corresponds to ^ 30 pc. 

2.2. AH03J1 725-34.4 

The second cluster, AH03 J1725-34.4, embedded in the star- 
forming region NGC 6357, is located at ^ 0?26 (or ^ 7.5 
pc in projection) to the southeast from Pismis 24 (Dias et al. 
12002b . The available information on this cluster is very scarce. 
AH03 J1725-34.4 contains at lea st four OB stars (^Neckel lT984l 
Damke, Barba & Morrell [2006b . Two of them, [N78] 49 and 
[N78] 51, are very massive ones, with spectral types (04 III((f*)) 
and 03.5 V, respectively; Damke et al. 12006b comparable to 
those of the most massive stars in Pismis 24. Unfortunately, 
AH03 J1725-34.4 was not covered by the Chandra observations 
of NGC 6357, so that the actual extent and the stellar content of 
the cluster are still unknown. The presence of two very massive 
stars in the cluster, however, suggests that it should be as young 
(~ 1 - 2 Myr) and at least half as massive as Pismis 24. 

3. Search for bow shocks 

A rich massive stellar content of Pismis 24 and the presence 
of very massive (binary) stars in both clusters associated with 
NGC 6357 suggest that this star-forming region was very ef- 
ficient in producing runaway stars (e.g. Gvaramadze 120071 
Gvaramadze, Gualandris & Portegies Zwart |2009al Gvaramadze 
& Gualandris 12011b . One can, therefore, expect to find some 
of the runaways via detection of their associated bow shocks - 
the natural attributes of supersonically moving stars (Baranov, 
Krasnobaev & Kulikovskii |1971[ Weaver et al. 11977b . The char- 
acteristic scale of a bow shock, /, depends on the number den- 
sity of the ambient medium, n, on the stellar mass-loss rate, 
M, and on the stellar space velocity, v*, as follows: / oc n~^^^, 
I oc M^/^ and / oc ^ (Baranov et al. 11971b . The bow shocks 
are usually most prominent in the mid-infrared (e.g. van Buren 
& McCray [T9881 van Buren, Noriega-Crespo & Dgani [T995t 



Paper I; Gvaramadze et al. I2010dl 1201 Ibb , but can also be de- 
tected in the optical (e.g. Kaper et al. |1997t Bro wn & Bomans 
2005; Paper I) and radio (Benaglia et al. 12010b wavebands. It 
should be noted that only a minority (^ 20 per cent) of run- 
away OB stars are associated with (detectable) bow shocks 
(van Buren [T9931 van Buren et al. [T9951 Hutholf & Kaper 
120021 Gvaramadze et al. l2010db . The paucity of the bow-shock- 
producing stars is mostly due to the fact that the majority of 
runaway stars are moving through a low-density, hot medium, 
so that the emission measure of their bow shocks is below the 
detection limit or the bow shocks cannot be formed at all be- 
cause the sound speed in the local interstellar medium is higher 
than the stellar space velocity (Kaper, Comeron & Barziv J999l 
Huthoff' & Kaper I2OO2I) . Moreover, the bow shocks generated by 
runaway stars receding from NGC 6357 and interacting with the 
dense material of the background (parent) molecular cloud (see 
Sect.[2l) would be very compact and therefore hardly detectable, 
while those projected against the Hii region might be hidden 
by its bright emission (see Fig.[T]). From this it follows that the 
actual number of stars ejected from the clusters in NGC 6357 
could several times exceed the number of detected bow- shock- 
producing stars. 

To search for bow shocks around NGC 6357, we used the 
archival data from the Mid-Infrared Galactic Plane Survey (Price 
et al. 2001), the 24 and 70 Micron Survey of the Inner Galactic 
Disk with MIPSEl (MIPSGAL; Carey et al. [2009b and the prelim- 
inary data release of the Mid-Infrared All Sky Survey carried out 
with the Wide-field Infrared Survey Explorer (WISE; Wright et 
al. 2010). The first survey, carried out with the Spatial Infrared 
Imaging Telescope onboard the Midcourse Space Experiment 
(MSX) satellite, covers the entire Galactic plane within \b\ < 5° 
and provides images at 1 8 arcsec resolution in four mid-infrared 
spectral bands centred at 8.3 yum (band A), 12.1 jim (band C), 
14.7 yum (band D), and 21.3 yum (band E). The MIPSGAL sur- 
vey (carried out with the Spitzer Space Telescope) mapped 278 
square degrees of the inner Galactic plane (\b\ < 1° is covered 
for 5° < / < 63° and 298° < / < 355° and \b\ < 3° is covered for 
|/| < 5°) and provides 24/dm images at 6 arcsec resolution. The 
current release of the WISE survey covers 57 per cent of the sky 
and provides images at four wavelengths: 3.4, 4.6, 12 and 22 //m, 
with angular resolution of 6.1, 6.4, 6.5 and 12.0 arcsec, respec- 
tively. The advantage of the MSX and WISE surveys is that they 
cover the whole Galactic plane and extend to higher Galactic 
latitudes than the MIPSGAL one, thereby allowing us to search 
for high- velocity runaways ejected at large angles to the Galactic 
plane (which could be off' the region covered by the MIPSGAL 
survey). On the other hand, the better angular resolution of the 
MIPSGAL survey could be vital for the detection of bow shocks 
generated by stars moving through the dense gas of the parent 
molecular cloud. To search for possible optical counterparts to 
the bow shocks and to identify their associated stars, we used 
the Digitized Sky Survey II (DSS-II; McLean et al. 2000). 

Assuming that the age of the star clusters in NGC 6357 is 
~ 1 - 2 Myr and adopting the typical velocity of runaway stars 
of ~ 30 - 50 kms"^ one finds that stars leaving the clusters at 
the very beginning of their dynamical evolution would be con- 
fined within -1-3° from NGC 6357. Thus, one can neglect the 
eff'ect of the Galactic gravitational potential on their trajectories. 
Correspondingly, one can expect that the bow shocks produced 
by the ejected stars would be directed away from their parent 



^ MIPS = Multiband Imaging Photometer for Spitzer (Rieke et al. 
[20041 . 
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Fig.l. 3° X 3° 21.3 Jim (MSX 
band E) image of the Hii 
region NGC 6357 (containing 
the star clusters Pismis 24 and 
AH03 J1725-34.4) and its en- 
vironments. The positions of 
seven bow shocks detected 
with MSX and Spitzer are 
marked by large and small 
circles, respectively. The po- 
sition of the two most mas- 
sive stars in Pismis 24 (which 
define the centre of the clus- 
ter) is indicated by a large 
diamond. The position of the 
most massive star in the cluster 
AH03J1725-34.4 is marked 
by a small diamond. The po- 
sition of IRAS 17207-3404 is 
indicated by a square, while 
the positions of three circular 
shells are indicated by crosses 
(see text for details). 



clusters, provided that there are no peculiar large-scale flows in 
the interstellar medium through which the stars are moving. 

First, we searched for bow shocks using the MSX data. 
The search was carried out in a 12° wide area elongated along 
the Galactic plane and centred at the longitude of NGC 6357 
(/ ^ 353°). Along the Galactic latitude, the search was lim- 
ited by the MSX coverage (\b\ < 5°), so that potentially we 
were able to detect bow shocks produced by stars leaving the 
cluster immediately after its formation and moving perpendicu- 
larly to the Galactic plane with a peculiar (tangential) velocity 
of ~ 100 km s"^. The visual inspection of MSX 21.3 jim images 
revealed five bow shocks (indicated in Fig. [T] by large circles). 
All of them have a clear arc-like structure that opens towards 
NGC 6357 (Figs. [2116]), which suggests that these structures are 
generated by stars expelled from this star-forming region. Bow 
shock 4 has a prominent gap on its leading edge (Fig.[?J see also 
below). All five bow shocks are most prominent at 21.3 yum, al- 
though some of them are visible in other MSX wavebands as 
well. Bow shock 3 has an obvious optical counterpart in the 
DSS-II (see Fig.|4]), while bow shocks 1 and 5 are embedded 
in H II regions (the H ii region associated with bow shock 5 is 
known as GAL 351.66+00.52; Lockman[T989]). 

Then, we used the data from the MIPSGAL survey. Four of 
the bow shocks discovered with MSX were covered by this sur- 
vey and we give their MIPS 24yum images in Fig.fT] We also 
discovered two new bow shocks whose orientation is consis- 





2'arcmin 



Fig. 2. Left: MSX 21 .3 yum image of bow shock 1 . Right: DSS-II 
(red band) image of the same field. The position of the bow- 
shock-producing star (star 1) is marked by a circle. 



tent with the possibility that their associated stars were ejected 
from NGC 6357 (see Figs.[8]and[9]and compare them with Fig.[Tl 
where the positions of these bow shocks are indicated by small 
circles). The angular extent of both bow shocks (~ 30 arcsec) 
is several times smaller than that of the bow shocks discovered 
with MSX, which can be understood if the stars generating these 
two bow shocks are either moving through the denser ambi- 
ent medium (e.g. the background molecular cloud) or/and have 
weaker winds and higher peculiar velocities. Neither of these 
bow shocks were detected in the optical range. 
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Fig. 3. Left: MSX2\3 yum image of bow shock 2. Right: DSS-II 
(red band) image of the same field. The position of the bow- 
shock-producing star (star 2) is marked by a circle. 
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Fig. 4. Left: MSX 21.3 yum image of bow shock 3. Right: DSS-II 
(red band) image of the same field. The position of the bow- 
shock-producing star (star 3 or HD 319881) is marked by a cir- 
cle. 
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Fig. 5. Left: MSX 8.3 yum image of bow shock4. Right: DSS-II 
(red band) image of the same field. The position of the bow- 
shock-producing star (star 4) is marked by a circle. 
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Fig. 6. L^/^.- M^X21.3 yum image of bow shock 5. Right: DSS-II 
(red band) image of the same field. The position of the bow- 
shock-producing star (star 5 or [N78] 34) is marked by a circle. 
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Fig. 7. From left to right, and from top to bottom: MIPS 24 yum 
images of bow shocks 1, 2, 3, and 5. The positions of the associ- 
ated stars are marked by circles. The orientation and the scale of 
the images are the same. 
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Fig. 8. Left: MIPS 24yum image of bow shock 6. Right: DSS-II 
(red band) image of the same field. The position of the bow- 
shock-producing star (star 6) is marked by a circle. 
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Fig. 9. Left: MIPS 24yum image of bow shock?. Right: DSS-II 
(red band) image of the same field. The position of the bow- 
shock-producing star (star 7) is marked by a circle. 



The use of the WISE data did not result in discovery of new 
bow shocks around NGC 6357. All seven bow shocks detected 
with MSX and Spitzer are clearly visible in the WISE data as 
well. In Fig.[TOlwe present images of bow shock 4 (not covered 
by the MIPS GAL survey) at all four WISE wavelenghts, show- 
ing its curious fine structure (most prominent at 12 and 22 yum). 
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Fig. 10. From left to right, and from top to bottom: WISE 22, 
12, 4.6 and 3.4 yum images of bow shock 4. The position of the 
associated star (star 4) is marked by a circle. The orientation and 
the scale of the images are the same. 
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Fig. 11. From left to right, and from top to bottom: IRAC 8, 4.5 
and 3.6 yum, and DSS-II (red band) images of the candidate bow- 
shock-producing star (star 8) embedded in the extended region 
of infrared emission (known as IRAS 17207-3404). The orien- 
tation and the scale of the images are the same. 



The possible origin of the gap in the leading edge of bow shock 4 
and the cirrus-like filaments around it are discussed in Sect.[6l 

A by-product of our search for bow shocks using the Spitzer 
data is the discovery of a compact arcuate nebula attached to 
one of the candidate OB stars revealed with Chandra (namely, 
star 9 from Table 7 of Wang et al. 2007). The MIPS 24yum im- 
age of this nebula is saturated, so that we give in Fig.[TT] the 
Spitzer 8, 4.5 and 3.6 yum images obtained with the Infrared 
Array Camera (IRAC; Fazio et al. 12004b within the Galactic 
Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE; 
Benjamin et al. 12003b . The orientation and the small angular 
size of the nebula (embedded in the more extended region of 
infrared emission known as IRAS 17207-3404; marked in Fig[T] 
by a square) are consistent with the interpretation of the neb- 
ula as a bow shock generated by a massive star ejected from 
Pismis 24 and moving through the dense background molecu- 
lar cloud. Correspondingly, the extended IRAS source can be in- 
terpreted as an H ii region created b y the ion izing emission of 
the star (cf. Paper I; Gvaramadze et al. l2011bl) . Alternatively, the 
nebula could be a circumstellar (toroidal) shell, similar to those 
observed around some evolved massive stars (see, e.g.. Fig. 4 in 
Gvaramadze, Kniazev & Fabrika l20 1 Oci f or the IRAC 5.8 yum im- 
age of a nebula around the candidate luminous blue variable star 
MN56). The 2MASS photometry of the star, however, suggests 
that it could be a main-sequence late 0-type star (see Sect.|4]), 
so that we consider the bow- shock interpretation of the arcuate 
nebula as more preferable. Note that the candidate OB star asso- 
ciated with the nebula is separated from the centre of Pismis 24 
by only ^ 9 arcmin, i.e. it is located (at least in projection) well 
within the cluster halo. For the sake of completeness, we give 
the details of this star (hereafter star 8) in Table[T] 

Another by-product of the search for bow shocks with 
Spitzer is the discovery of three ring-like shells with central 
point sources. We present these shells and discuss their origin 
in Sect.[5l 



4. Bow-shock-producing stars 

Inspection of the DSS-II (red band) images suggested that the 
bow shocks discovered with MSX and Spitzer are associated with 
stars indicated in Figs. [2] -[9] by circles. The details of these stars 
are given in Table[Tl Columns 2, 3, 5, and 6 give the equato- 
rial coordinates of the stars and their /- and -magnitudes (all 
taken from the 2MASS catalogue; Skrutskie et al. 120061) . while 
column 4 gives the visual magnitudes of the stars taken from 
the NOMAD catalogue (Zacharias et al. 2004). Columns 7 and 
8 give the /sTs-band interstellar extinction towards the stars and 
their absolute magnitudes in the same band derived under the 
assumption that the stars are located at the same distance as 
NGC 6357. Column 9 provides spectroscopic and photometric 
(enclosed in brackets) spectral types of the stars (see below). 

Using the SIMBAD databaseQ and the VizieR catalogue ac- 
cess toojj, we found that two of the five bow shocks discov- 
ered with MSX are generated by known OB stars (with spec- 
troscopically derived spectral types). Namely, it turns out that 
bow shocks 3 and 5 are produced by the 06 Vn (Drilling & 
PerryUMB star HD 319881 and the 08/BOIII/V (Neckel[M3 
star [N78] 34, respectively. To determine the nature of the other 
three stars (hereafter stars 1, 2 and 4), we obtained their spectra 
with the South African Astronomical Observatory (SAAO) 1.9- 
m telescope (see Sect. l4.1l) . The stars associated with the bow 
shocks discovered with Spitzer (hereafter stars 6 and 7) still have 
to be observed spectroscopically. If, however, we assume that 
these two stars are located at the same distance as NGC 6357, 
then one can estimate their spectral types using the 2MASS 
photometry and calibration of absolute magnitudes and intrin- 
sic colours (see Sect. l4.2l) . 



^ http : //simbad. u- strasbg . fr/simbad/ 

^ http : //web viz . u- strasbg . f r/viz-bin/VizieR 
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Table 1. Details of seven bow- shock-producing stars (stars 1-7) and one candidate bow- shock-producing star (star 8) detected 
around NGC 6357. The last column provides spectroscopic and photometric (enclosed in brackets) spectral types of the stars (see 
text for details). 



Star 


RA (2000) 


Dec. (2000) 


V 


/ 








Spectral type 


1 


17 27 11.23 


-34 14 34.9 


11.77 


8.42 


7.84 


0.52 


-3.84 


07.5:^ (07 V) 


2 


17 22 03.43 


-34 14 24.1 


13.59 


8.94 


8.01 


0.75 


-3.89 


05.5:^ (06.5-7 V) 


3 or HD 319881 


17 28 21.67 


-34 32 30.3 


10.14 


7.54 


7.10 


0.43 


-4.48 


06 Vn^ (05 V) 


4 


17 18 15.40 


-34 00 06.1 


11.80 


8.19 


7.53 


0.57 


-4.19 


06.5:^ (06 V) 


5 or [N78] 34 


17 22 05.62 


-35 39 55.5 


13.11 


8.68 


7.82 


0.71 


-4.04 


O8/B0IIW^ (06.5 V) 


6 


17 22 50.02 


-34 03 22.4 


13.90 


9.68 


8.81 


0.71 


-3.05 


(BOV) 


7 


17 27 12.53 


-33 30 40.0 


12.18 


9.21 


8.65 


0.51 


-3.01 


(BOV) 


8 


17 24 05.62 


-34 07 09.5 


11.74 


9.03 


8.46 


0.51 


-3.20 


(09-9.5 V) 



Notes. this work. Drilling & Perry ([19811 ). Neckel ([T9841) . 



4.1. Spectral types of the bow-shock-producing stars 1, 2 
and 4 

The observations of stars 1 , 2 and 4 were performed with the 
Cassegrain spectrograph using a long slit of 3'xl.8'' on 2009 
June 1 and 2. Grating with 300 lines mm"^ was used with the 
266x1798 pixels CCD detector SITe with two setups with the 
wavelength coverage of ~ 3760 - 7750 A for stars 1 and 4 and 
~ 3560 - 7550 A for star 2. For both setups the spectral reso- 
lution was ~2.3 A pixel" ^ (FWHM ~7 A) along the dispersion 
direction and the scale along the slit was Vf4 pixel" ^. The see- 
ing during the observations varied from 1.0 to 2.1 arcsec from 
object to object, but was stable for each object. Exposure times 
were 3 x 300, 3 x 600 and 3 x 300 sec for stars 1, 2 and 4, respec- 
tively. The observations of star 2 were made in variable condi- 
tions, with drastic drops in transparency, resulting in the worse 
quality of the resulting spectrum. Spectra of Cu-Ar comparison 
arcs were obtained to calibrate the wavelength scale. The spec- 
trophotometric standard stars LTT7379 and LTT 3864 were ob- 
served with a slit width of 5 arcsec at the beginning and at the 
end of the nights for flux calibration. Data reduction was per- 
formed using standard procedures (see, e.g., Kniazev et al. l2QQ8l 
for details). 

Figure[T2]presents the resulting normalized spectra of stars 1 , 
2 and 4 in the AA4000 - 7000 A region. All three spectra show 
numerous narrow absorption lines of hydrogen. He i and He ii, 
which is typical of 0-type stars. Siiv /14089 and Oiii /15592 ab- 
sorptions are also visible in the spectra. We detected the emis- 
sion of C III /15696 in all three spectra (which characterizes the 
stars of the 04-08 spectral type; Walborn 11980b , although in 
the spectrum of star 2 the detection is marginal. No forbidden 
lines can be seen in the spectra. Numerous difl'use interstellar 
bands (DIBs) are present in the spectra. A blend of strong Nai 
/15890, 96 absorption lines is of circumstellar or interstellar ori- 
gin, while the strong absorptions visible at A > 6800 A are tel- 
luric in origin. 

To determine the spectral types of stars 1 and 4, we used 
the traditional classification criteria based on the ratio of the 
equivalent widths (EWs) of Hi A4 411 and Hen /1454 1 line s 
(Morgan, Keenan & Kellman [19431 C onti & Alschuler [19711). 
Using Table 3 in Conti & Alschuler (11971b and EWs given 
in Table[2l we found the spectral types of stars 1 and 4 to be 
07.5 and 06.5, respectively, with an uncertainty of ±1 subtype. 
The more earlier subtype inferred for star 4 is consistent with 
the increase of the Cm /15696 emission line strength with lu- 
minosity ( Walbo rn 11980b . Similarly, using Table 5 in Conti & 
Alschuler iWUh and EWs of Siiv .14089 and Hei .14143 lines 
from Table[2l we found that both stars are supergiants. This in- 



ference, however, is inconclusive because of the large errors in 
the measurements of the Si iv /14089 line, which allow the two 
stars to belong to any luminosity class (see also below). 

The poor observational conditions and the relatively high ex- 
tinction towards star 2 (see Table[T]) degraded the spectrum of the 
star, which precludes us from using traditional classification cri- 
teria. Instead, we utilized the classification scheme for O stars 
in the yellow-green (/1/14800 - 5420 A) proposed by Kerton, 
Ballantyne & Martin (11999b , which is based on the ratio of EWs 
of the lines He i .14922 and He ii .1541 1 : 

SpT = (9.04 ± 0.10) -h (4.10 ± 0.23) log EW4922/EW5411 . (3) 

Using Eq. ([3]) and EWs given in Tabled we found SpT^ 5.5. 
Similarly, applying this equation for stars 1 and 4 gives their 
spectral types of 08 and 07, respectively, which agree well with 
the estimates based on the traditional classification criteria. The 
spectral type can also be estimated using separately EW4922 and 
EW5411 (Kerton et al l 19991) : 

SpT = (4.82 ± 0.03) (7.92 ± 0.56)EW4922 , (4) 

SpT = (12.80 ± 0.41) - (7.08 ± 0.58)EW54ii . (5) 

From Eqs. (O and ([51), we found SpT^ 6.0 and 5.0, respectively; 
both estimates are consistent with that derived from Eq. ([3]). Still, 
although it is clear that star 2 is of 0-type, the estimates of its 
spectral type suff'er from the large errors of measured EWs (see 
Table[2l). 

The lack of an exact spectral classification does not allow 
us to compare the spectroscopic distances to stars 1 , 2 and 4 
with the distance to NGC 6357. Instead, we use the assumption 
that these stars were ejected from NGC 6357 to constrain their 
spectral types and luminosity classes. Using Eq. ([T]) and adopt- 
ing the intrinsic colou r, for O stars, (/ - K^)o = -0.21 mag 
(Martins & Plez 12006b , we calculated the ^s-band extinction, 
Ak^ , towards these stars and their absolute magnitudes Mk^ (see 
Table[T]). Then using the calibration of absolute infrared magni- 
tudes for Galactic O stars by Martins & Plez (2006), we found 
that all three stars should be on the main- sequence with spectral 
types of 07 V, 06.5-7 V, and 06 V, which agree well with those 
derived from spectroscopy (see Table[Tl where the photometric 
spectral types are enclosed in brackets). 

4.2. Constraining the spectral types of the other 
bow-shock-producing stars 

Similarly, assuming that the other bow- shock-producing stars 
around NGC 6357 are located at the same distance as this Hii 
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Fig. 12. Normalized spectra of stars associated with bow shocks 4 (top), 1 (middle) and 2 (bottom) with principal lines and most 
prominent DIBs indicated. The spectra of stars 4 and 2 are shifted upwards and downwards by ^ 0.3 continuum flux unit, respec- 
tively. 

Table 2. Equivalent widths (EWs; in A) of absorption lines used for spectral classification of bow-shock-producing stars 1, 2 and 4. 



Star 


Si IV A40S9 


He 1^4143 


He 1/14471 


Hen /14541 


He 1^4922 


Hen i5411 


1 


0.62+0.35 


0.30+0.15 


0.90+0.12 


0.75+0.09 


0.58+0.08 


1.07+0.05 


2 










0.15+0.14 


1.07+0.25 


4 


0.53+0.34 


0.14+0.06 


0.70+0.12 


1.00+0.12 


0.34+0.13 


1.15+0.05 



region, we calculated Ak^ and Mk^ for these stars and estimated 
their photometric spectral types (see Table[T]). For stars 3 and 5 
we found spectral types of 05 V and 06.5 V, respectively. The 
former estimate agrees reasonably well with the spectroscopi- 
cally derived spectral type, while the latter one is inconsistent 
with the spectral classification reported by Neckel (119841) . We 
note, however, that the spectral types given by Neckel (119841 ) 
are systematically later than those derived from th e mor e re- 
cent observations (Massey et al. 120011 Damke et al. 12006b . For 
example, one of the brightest (and most massive) stars in the 
cluster AH03J1725-34.4, [N78]51, was classified by Neckel 
(119841 ) as an 09 star, while in the recent study of this cluster by 
Damke et al. (120061) the spectral type of the star was found to be 
03.5 V((f)). 

Applying the same procedure for stars 6 and 7 suggests that 
both of them could be of the same spectral type of BO V (so 
that the smaller size of the bow shocks generated by stars 6 and 
7 could be because of weaker winds of these stars). And fi- 



nally, we estimated the spectral type of the candidate bow- shock- 
producing star (star 8) to be 09-9.5 V. Additional spectroscopic 
and photometric observations of the bow- shock-producing stars 
around NGC 6357 are warranted. 



4.3. Proper motions and likely parent clusters of 
bow-shock-producing stars 

In Table[3] we provide two representative proper motion mea- 
surements for each bow- shock-producing star, found with 
VizieR. For most stars we used the measurements from the 
most recent catalogues, namely UCAC3 (Zacharias et al. 12010b 
and PPMXL (Roser, Demleitner & Schilbach [20T0b . When the 
UCAC3 proper motion is marked "with doubts" (because it re- 
lies on "less than two good matches"), we used the data from the 
NOMAD catalogue (Zacharias et al. 2004 ). Although most mea- 
surements are insignificant (i.e. the measurement uncertainties 
are comparable to the measurements themselves), it is of inter- 
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Table 3. Proper-motion measurements for the seven bow- shock-producing stars (stars 1-7) and one candidate bow- shock-producing 
star (star 8) around NGC 6357. Two measurements are given for each star to indicate the uncertainties in the measurements. For 
each data set, the peculiar (transverse) velocities (in Galactic coordinates) were calculated and added to the table. A likely parent 
cluster for each star is indicated in the last column. 



Star 


jJLa COS 6 




Ref. 


Vl 




Parent cluster 




mas yr~^ 


mas yr~^ 




kms~^ 


kms~^ 




1 


-1.5 + 1.9 


-1.6 + 2.0 


1 


-4.4 + 16.1 


9.9 + 15.3 


AH03J1725-34.4 


1 


-3.8 + 3.0 


-5.0 + 3.0 


2 


-7.1 +24.2 


10.0 + 24.2 




2 


-8.5 + 4.7 


-3.2 + 4.7 


3 


-46.8 + 37.9 


48.5 + 37.9 


Pismis 24 


2 


-8.9 + 12.2 


-11.3 + 12.2 


2 


-102.4 + 98.3 


14.1+98.3 




3 (HD 319881) 


0.0+ 1.4 


-2.8 + 1.0 


1 


-5.7 + 8.9 


-5.3 + 10.5 


Pismis 24 


3 (HD 319881) 


3.4 + 2.6 


-3.2 + 2.7 


2 


6.8 + 21.8 


-29.9 + 21.0 




4 


-4.4 + 2.0 


0.9 + 1.1 


1 


-1.0+11.8 


39.7 + 10.4 


Pismis 24/AH03 J1725-34.4 


4 


-9.0 + 2.9 


2.3 + 3.0 


2 


-13.2 + 23.9 


76.5 + 23.6 




5 ([N78] 34) 


-4.9 + 4.8 


11.7 + 4.8 


3 


69.0 + 38.7 


93.1 +38.7 


AH03J1725-34.4 


5 ([N78] 34) 


-11.8 + 12.2 


18.7 + 12.2 


2 


83.9 + 98.3 


170.9 + 98.3 




6 


-4.6 + 3.8 


0.8 + 2.1 


1 


-2.5 + 22.3 


41.0 + 27.0 


AH03J1725-34.4 


6 


-8.1 + 12.2 


0.1 + 12.2 


2 


-23.1 +98.3 


61.1+98.3 




7 


-3.9 + 2.1 


0.0+1.9 


1 


-4.7 + 15.8 


33.1 + 16.4 


Pismis 24 


7 


-7.3 + 3.0 


-3.2 + 3.0 


2 


-41.5 + 24.2 


41.5 + 24.2 




8 


-6.0+1.8 


1.9 + 2.6 


1 


-1.4+19.1 


55.5 + 16.8 


Pismis 24 


8 


-9.4 + 3.0 


-2.5 + 3.0 


2 


-46.2 + 24.2 


58.1 +24.2 





References. (1) UCAC3 (Zacharias et al. [20T0l) : (2) PPMXL (Roser et al. [20T0l ): (3) NOMAD (Zacharias et al. 12004] ). 



est to check whether or not the orientation of the stellar peculiar 
velocities implied by these measurements is consistent with the 
distribution of the stars on the sky and with the orientation of 
their associated bow shocks. 

To convert the observed proper motions into the transverse 
peculiar velocities of the stars, we used the Galactic constants 
Rq = 8.0 kpc and ©o = 240 km s'^ (Reid et al. 2009) and 
the solar peculiar mo tion (U q, V©, Wq) = (10.0, 11.0, 7.2) kms"^ 
(McMillan & Binnev l2QTQl) . The derived velocity components in 
Galactic coordinates are given in columns 5 and 6 of Tabled For 
the error calculation, only the errors of the proper motion mea- 
surements were considered. Obviously, the orientation of pecu- 
liar velocities of stars 1, 2, 3, 4, 6 and 8 are consistent (within 
their margins of error) with the orientation of their bow shocks 
and with the possibility that these stars are moving away from 
NGC 6357. For stars 5 and 7 we found that their peculiar ve- 
locities have "wrong" orientations, i.e. are inconsistent with the 
orientation of the associated bow shocks. 

Because of the low significance of proper motion measure- 
ments, we rely only on the orientation of the bow shocks in the 
identification of parent clusters for bow- shock-producing stars. 
We caution, however, that the close proximity of the clusters on 
the sky and the absence of proper motion measurements for the 
clusters make the identification ambiguous. Thus, the likely par- 
ent clusters indicated in the last column of Table[3] should be 
considered as tentative. 

5. Three circular shells 

As we mentioned in Sect.O one of the by-products of the search 
for bow shocks around NGC 6357 with Spitzer is the discov- 
ery of three circular shells. The positions of these shells on the 
sky are marked in Fig.[T]by crosses. The appearance of all three 
shells (see Figs.[T3l-[T5]) is typical of circumstellar shells created 
by luminous blue variable (LBV) stars and Wolf-Rayet stars of 
late WN-types (WNL) (see Gvaramadze et al. 1201 Pel for Spitzer 
images of such shells). 

One of the shells (hereafter shell 1) is located at only ^ 0?6 
(or ^ 18 pc in projection) to the northwest of NGC 6357. 




Fig. 14. Left: MIPS 24 yum image of shell 2. Right. IRAC 8 yum 
image of the same field. 



The MIPS 24yum image of this shell (see left panel of Fig.[T3]) 
shows a clear ring-like structure with a diameter of ^ 25 arcsec 
and enhanced brightness along the southeast rim faced towards 
NGC 6357. Fig. [O] also shows the central point-like source be- 
ing smeared-out somewhat because of the 6 arcsec resolution of 
the MIPS 24yum image. We interpret this source as a massive 
evolved star physically related to the shell. Many dozens of sim- 
ilar evolved stars were recently revealed through the detection 
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Fig. 15. Left and Middle: MIPS 24 jim and WISE 22 jim images of shell 3, respectively. The MIPS image is truncated because the 
MIPSGAL survey covers only a part of the shell. Right. DSS-II (red band) image of the same field. The central star of shell 3 is 
marked by a circle. 



Table 4. Circular shells around NGC 6357 and their central stars. 



Shell 


RA (2000) 


Dec. (2000) 


V 


[3.6] 


[8.0] 


Diameter (arcsec) 


1 


17 22 36.7 


-33 49 10 








25 


2 


17 26 19.11^ 


-35 32 37.7^ 




13.67^ 


8.29^ 


40 


3 


17 20 31.76^ 


-33 09 48.9^ 


143^ 






210 



Notes. ^""^ GLIMPSE Source Catalog (lhttp://irsa.ipac?caltech.edu/applications/Gator/ ). NOMAD (Zacharias et al. l2004l) . 



of their infrared circumstellar shells with Spitzer and follow-up 
spectroscopy (Gvaramadze et al. |2009bl I2010al bx: Wachter et 
al. 2010. 2011j). The DSS-II and 2MASS images, however, do 
not show any counterpart to the putative central star of shell 1 . 
The most plausible explanation of the absence of the 2MASS 
counterpart is that the central star is highly absorbed by the fore- 
ground dusty material. In this connection, we note that the cen- 
tral stars of some circumstellar shells presented in Gvaramadze 
et al. (1201 Oct also cannot be seen in the 2MASS images, while 
they are visible in the IRAC (3.6, 4.5, 5.8 and 8.0 jim) ones. 
Unfortunately, shell 1 was not covered by the GLIMPSE sur- 
vey, so that at the moment we cannot prove the existence of the 
central source that is apparently visible at 24 jum. Deep infrared 
imaging of shell 1 is therefore highly desirable. Detection of the 
central star and its follow-up spectroscopy would allow us to un- 
veil the nature of the shell and potentially link the star to the 
star-forming region NGC 6357. 

If one assumes that shell 1 is produced by a massive star 
ejected from NGC 6357, then the linear size of the shell is ^ 0.2 
pc, which is comparable to that of some LBV shells (e.g. Buemi 
et al. 12010b . If, however, the actual size of the shell is 1-2 pc 
(which is more typical of circumstellar shells associated with 
LBV and WNL stars), then the distance to this shell is ^ 7 - 14 
kpc. At this distance the star would be at a height of ^ 180 - 360 
pc above the Galactic plane, which would imply that it is a run- 
away. 

Eor the limiting magnitude of the 2MASS survey of 14.8 
and the K^-hsind absolute magnitude of LBV and WNL s tars of 
^ -9 and -6 m ag, res pectively (Blum, Depoy & Sell gren 119951 
Crowther et al. 12006b , one finds that Ak^ ^ 10-13 mag (or 
Ay ^ 90-1 10 mag) if the star is located at 1.7 kpc. An extinction 
this high could be understood if the star were located within (or 
behind) the dense molecular cloud associated with NGC 6357 
(see Sect.[2l). 



The second shell (hereafter shell 2) is located at ^ 1?3 to 
the south of NGC 6357. Eig.IH shows the MIPS 24yum (right 
panel) and IRAC 8 yum images of this shell. The ring-like struc- 
ture of shell 2 is more obvious in the 8 jim image, from which 
we measured a diameter of the shell of ^ 40 arcsec (^ 0.3 pc 
at a distance of 1.7 kpc). The central star of shell 2 is visible not 
only at 24yum, but also in all four IRAC bands. We also found 
a weak counterpart to this star in the 2MASS K^-hmid image, 
which, however, is not Hsted i n the 2 MASS All-Sky Catalog of 
Point Sources (Skrutskie et al. 120061) . 

The third shell (hereafter shell 3) is located at ^ 1?4 to the 
northwest of NGC 6357. Unfortunately, the MIPSGAL survey 
covers only a part of this quite extended (^ 3.5 arcmin in di- 
ameter) shell with a prominent blow-up to the northeast (see 
left panel of Eig.[T5]). The full extent of shell 3, however, is vis- 
ible in the WISE 22yum image (middle panel of Fig. [15]), which 
shows an almost circular limb-brightened nebula and its central 
star (visible also in the DSS-II image; right panel of Fig.[T5l). 
Using the SIMBAD database, we identified the central star with 
the emissio n-line star Hen 3-1383, classified in the literature as 
an Me (The [T96T[ Sanduleak & Stephenson [1971 Allen [19781) 
or a candidate symbiotic (Henize [1976b star. Our recent optical 
spectroscopy of Hen 3-1383 with the Southern African Large 
Telescope (SALT), however, showed that the spectrum of this 
star is almost identical to those of the prototype LBV P Cygni 
(see, e.g., Stahl et al. 11993b and the recently discovered candi- 
date LBV MN 112 (Gvaramadze et al. l2QlQbb , which impHes the 
LBV classification for Hen 3-1383 as well (Gvaramadze et al., in 
preparation). 

The details of the three shells are summarized in Tabled In 
this table we give approximate coordinates of the geometric cen- 
tre of shell 1 and the coordinates of the central stars associated 
with two other shells. For the central star of shell 2 we give its 
IRAC 3.6 and 8.0 jim magnitudes (taken from the GLIMPSE 
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Source Catalog, while for the central star of shell 3 we give its 
visual magnitude from the NOMAD catalogue. The last column 
of the table contains the angular size of the shells. 

6. Discussion 

We explored the idea that young star clusters lose a significant 
fraction of their massive stars because of gravitational few-body 
interactions between the cluster members. The efficiency of dy- 
namical star ejection depends mainly on the star number den- 
sity in the cluster core, which is maximum either already during 
cluster formation or becomes maximum later on owing to the 
Spitzer's mass segregation instability (Spitzer [T969l) . Mass seg- 
regation could be a rapid 0.5 Myr) process if the cluster was 
born in a compact configuration (e.g. Kroupa lTOOSi) , so that dy- 
namical star ejections could be effective well before the cluster 
starts to expel its members via binary- supernova explosions (the 
latter process becomes efficient several Myr after the formation 
of the first massive binary stars). 

The runaway OB stars spread all around the parent cluster 
constitute the population of massive field stars. It is believed 
that most, if not all, field OB stars were formed in the clustered 
way (Lada & Lada [20031 Smith, Longmore & Bonnell ^0^ 
and subsequently found themselves in the field either because 
of the dynamical processes in the parent clusters or because of 
rapid cluster dissolution (e.g., Kroupa & Boilv [20021 de Wit 
et al. [20051 Schilbach & Roser [20081 Gvaramadze & Bomans 
I2008bt Pflamm-Altenburg & Kroupa [20T0l) . Detection of mas- 
sive stars in the field around young clusters and linking them to 
these clusters is therefore important not only for the problem of 
the origin of field OB stars, but also for the problem of massive 
star-formation in general. 

The onset of OB star ejection from a young cluster very 
much depends on whether the massive stars form as tight bina- 
ries near the cluster centre or not (Clarke & Pringle[T992j). Thus 
N-body experiments can be performed with different assump- 
tions on the OB stellar population. To see which of the theo- 
retical (numerical) models of cluster dynamical evolution better 
correspond to observations, one should identify as many stars 
ejected from a given star cluster as possible. The detected run- 
away stars would constitute a sample for the future high-precise 
proper motion and parallax measurements with the space as- 
trometry mission Gaia, which are required to determine the tim- 
ing of ejections (e.g. Paper I) and thereby to constrain the initial 
conditions for young star clusters (Kroupa 2008) and to distin- 
guish between the primordial and the dynamical origin of mass 
segregation in these clusters. 

In the absence of reliable proper motion measurements for 
distant and/or highly obscured field OB stars, the only viable 
possibility to prove their runaway status and to identify their par- 
ent clusters is through the detection of bow shocks that are pos- 
sibly generated by these starfl The characteristic scale of bow 
shocks is determined (for the given number density of the ambi- 
ent interstellar medium and the peculiar velocity of the star) by 
the strength of the stellar wind, which makes the massive stars 
prone for producing detectable bow shocks. On the other hand, 
the fact that dynamical star ejection is most efl&cient in mass- 

^ Available at "http ://irsa.ipac .caltech.edu/applications/Gator/ 
^ The runaway status of some field OB stars can also be confirmed via 
measurement of their high peculiar radial velocities, but this channel for 
detection of runaways does not allow us to determine the star's direction 
of motion on the sky and therefore cannot be used to identify their birth 
clusters. 



segregated clusters (Oh et al., in preparation) implies that during 
the early dynamical evolution of these clusters the ejected stars 
are preferentially the most massive ones. This inference is con- 
sistent with the observational fact that the percentage of O stars 
among the runaways is higher than that of B stars (e.g. Stone 
119911 ) and could also be derived from the existence of an ex- 
tended (tens of parsecs) halo of very massive (02-type) stars 
around the very young (~ 1-2 Myr) star cluster R136 in the 
Large Magellanic Cloud (e.g. Brandl et al. 120071 Walborn et al. 
[2003) . 

Our search for OB stars running away from the star-forming 
region NGC 6357 has resulted in the discovery of seven bow 
shocks. The morphology of the bow shocks suggests that 
the associated stars were ejected from clusters Pismis24 and 
AH03 J1725-34.4, which are embedded in NGC 6357. The ac- 
tual number of ejected stars could be much higher because only 
a small fraction of runaway stars produce observable bow shocks 
(see Sect.[3l). Four of five (or even all; see Sect. l4.2l) bow-shock- 
producing stars with spectroscopically determined spectral types 
turn out to be 0-type stars, which provides additional evidence 
that young massive star clusters expel preferentially the most 
massive stars. 

Follow-up spectroscopy of three bow- shock-producing stars 
detected around NGC 6357 showed that all of them are 0-type 
stars. Thus, the search for bow shocks around young clusters 
and follow-up spectroscopy of their associated stars is a useful 
tool for detection of highly obscured, yet unknown OB stars. 
In some cases, the bow- shock-producing stars are not necessar- 
ily high- velocity runaways. For example, high- velocity outffows 
driven from young star clusters by the collective effect of stellar 
winds and ionizing emission can create bow shocks around (low- 
velocity; not necessarily massive) stars in the cluster's halo. In 
this case, the bow shocks are facing towards the cluster's centre. 
Numerous examples of such bow shocks were detected around 
young clusters and OB associations: Orion Cluster (Bally, ODell 
& McCaughrean[2000|, Trumpler 14 (Ascenso et al. [2Qn7b . M17 
and RCW 4 9 (Pov ich et al. [2Qn8]) . Cyg 0B2 (Kobulnicky, Gilbert 
& Kiminki [2010b . One can also imagine the situation when a 
slowly (~ 10 kms"^) moving massive star encounters a dense 
medium (e.g. the remainder of the parent molecular cloud) and 
its wind starts to interact with the ionized gas outflowing from 
the cloud surface towards the star with a velocity comparable 
to the sound speed (^ 10 kms"^). In the reference frame of 
this gas, the stellar motion is supersonic so that a bow shock 
is formed ahead of the (slowly) moving star. We speculate that 
bow shock 1 was probably generated in this way. This possibil- 
ity is suggested by the small separation (^ 10 pc in projection) 
of starl from its likely birth cluster AH03 J1725-34.4, which 
implies that the peculiar velocity of the star could be as low as 
5-10 km s"\ provided that it was ejected from the cluster -1-2 
Myr ago (cf. Table[3l). Moreover, the right panel of Fig. [2] shows 
that bow shock 1 and its associated star are located at the periph- 
ery of the H II region, from which one can infer that star 1 has 
met and ionized a density enhancement on its way. 

Finally, we discuss the origin of the gap at the leading edge 
of bow shock 4 and the curious cirrus-like filaments around it. In 
Gvaramadze et al. ( 201 lb), we suggested that similar cirrus-like 
filaments around the bow shock generated by the hi gh-mas s X- 
ray binary VelaX-1 (see Fig. 2 in Gvaramadze et al. l2011bl ) are 
caused by interstellar dust grains aligned with the local interstel- 
lar magnetic field and heated by the radiation of VelaX-1. This 
suggestion implies that bow shock 4 propagates along the local 
interstellar magnetic field, which affects the compression ratio of 
the shocked gas by making it less dense along the flanks of the 
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Fig. 16. MIPS 24 fim image of the bow shock generated by the 
BC0.7 la star HD 2905 and cirrus-like filaments around it (see 
text for details). 

bow shock (e.g. Draine & McKee [1993b . The density asymme- 
try in turn could aff'ect the heating of the dust grains accumulated 
behind the front of bow shock 4 and thereby could be responsi- 
ble for its brightness asymmetry. In the absence of reliable data 
on parameters of bow shock 4 and its associated star we refrain 
from a more detailed discussion of the problem. We note, how- 
ever, that similar gaps are typical of many other bow shocks. 
One of them (although less impressive and not at the symmetry 
axis of the bow shock) is shown in Fig. [161 which presents the 
MIPS 24 Jim image (Program Id.: 30088, PI: A.Noriega-Crespo) 
of a bow shock produced by the ECO. 7 la (Walborn 1972) star 
HD 2905. Like bow shock 4, the bow shock around HD 2905 is 
surrounded by cirrus-like filaments, some of which are appar- 
ently attached to the front of the bow shock near the gap. To 
conclude, if our suggestion on the origin of gaps is correct, then 
the cirrus-like filaments around bow shocks could serve as trac- 
ers of the local magnetic field, while the bow- shock-producing 
stars could serve as probes of the Galactic magnetic field. 

7. Summary 

Search for bow- shock-producing stars around the star-forming 
region NGC 6357 with two embedded young (~ 1 - 2 Myr) 
star clusters Pismis 24 and AH03 J1725-34.4 resulted in the dis- 
covery of seven bow shocks. The geometry of the bow shocks 
strongly suggests that all seven associated stars were ejected 
from NGC 6357. Because only a small fraction of runaway stars 
produce detectable bow shocks, the actual number of massive 
stars ejected from NGC 6357 could several times exceed the 
number of detected bow- shock-producing stars. This inference is 
consistent with the theoretical expectation that young star clus- 
ters lose a significant fraction of their massive star content at the 
very beginning of their dynamical evolution. Follow-up spec- 
troscopy of three detected bow- shock-producing stars showed 
that all of them are 0-type stars. This result is consistent with the 



observational fact that O stars prevail among the runaways and 
implies that search for bow shocks around massive clusters pro- 
vides a useful tool for discovery of new OB stars. A by-product 
of our search for bow shocks is the detection of three circular 
shells around NGC 6357, whose morphology is typical of LBV 
and WNL stars. Follow-up spectroscopy of the central star of 
one of these shells showed that its spectrum is almost identical 
to that of the prototype LBV P Cygni, which implies the LBV 
classification for this star as well. 
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